subroutine comput_flux2
use var
implicit none

    call comput_Qy_m
    call comput_Qy_p
    call L_F_local_y
    
end subroutine